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ABSTRACT 

Electric charge neutrality in a compact star provides an important relation- 
ship between the chiral dynamics and neutrino propagation in the star. Since the 
sudden drop of the electron density at the critical point of the first-order chiral 
phase transition, the oscillation for low energy neutrinos is significant and can 
be regarded as a signature of chiral symmetry restoration in the core of the star. 

Subject headings: dense matter — neutrinos — stars:neutron 

1. Introduction 

It is generally believed that there are two important QCD phase transitions in hot and 
dense nuclear matter. One of them is related to the deconfinement process in moving from 
a hadron gas to a quark-gluon plasma, and the other one describes the transition from the 
chiral symmetry breaking phase to the phase in which it is restored. Furthermore, it ap- 
pears from lattice simulations of QCD that both transitions coincide at a temperature of 
about T = 165 MeV at zero baryon density (Hwa 1990). The transitions may happen in the 
early universe and the early stage of relativistic heavy-ion collisions where the temperature 
is expected to be extremely high. The other often considered systems to realize the transi- 
tion condition are compact stars where the temperature is low but the baryon density can 
reach several times the normal nuclear density (Glendenning 2000). From the recent lattice 
QCD simulations, the phase structure at high density is much more rich than that at high 
temperature. Apart from the quark deconfinement and chiral restoration, there may be pion 
superfluidity (Kogut & Sinclair 2002;He,Jin,& Zhuang 2005) and even color superconduc- 
tivity (Alford,Rajagopal,& Wilczek 1998)in compact stars. Different from the temperature 
effect which leads to a second order phase transition or even only a crossover, the QCD phase 
transitions at high baryon density are all of first order (Karsch 2002) which can help us to 
extract signatures in experiments easily. 
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The signatures of a QCD phase transition that happened in a compact star depends 
strongly on the structure of the star, namely the equation of state of the star. Normally, a 
compact star is considered to be completely in a new state of matter, and a bag constant is 
needed to balance the pressure on the surface of the star (Glendenning 2000). In this case, 
there is no phase transition inside the star, and the signatures reflect only the properties of 
the new state of matter. If a phase transition occurs inside a compact star, the sudden change 
in the equation of state at the critical point of the first order phase transition may lead to 
some easily observable signatures. These signatures are connected to the phase transition 
itself and different from those related only to the new state of matter. 

Neutrinos play an important role in the study on compact stars (Bethe 1990;, Qian et 
al.l993;Prakash et al.l997;Yakovlev et al.2001). In this paper, we investigate the neutrino 
conversion in a compact star in the frame of an effective QCD model with chiral symmetry. 
While it is generally believed that the neutrino oscillation is difficult to occur in compact 
stars (Bethe 1990), the situation may be changed when the transition from chiral restoration 
phase to chiral breaking phase is taken into account. We will show that a sudden drop of the 
electron density associated with the first order chiral phase transition leads to an remarkable 
neutrino oscillation. 

The paper is organized as follows. In Section [2] we describe the chiral thermodynamics in 
the Nambu-Jona-Lasinio (NJL) model in mean field approximation and then by considering 
the electric charge neutrality constraint we associate the chiral properties with the inner 
structure of a compact star via the Tolman-Oppenheimer-Volkoff (TOV) equation. In Section 
[3] we study the effect of chiral phase transition on the neutrino oscillation through the 
neutrino propagation in the star. We summarize and discuss the obtained results in Section 
HI We use the natural unit of c = H = kg = 1 through the paper. 

2. Chiral Thermodynamics 

The study on QCD phase diagrams depends on lattice QCD calculations and effective 
models with QCD symmetries. Since there is not yet precise lattice result at finite baryon 
density due to the fermion sign problem(Karsch 2002), the structure of chiral phase transition 
at moderate baryon density is mainly investigated in many low energy effective models. 
One of the models that enables us to see directly how the dynamical mechanisms of chiral 
symmetry breaking and restoration operate is the NJL model(Nambu & Jona-Lasinio 1961) 
applied to quarks (Vogl & Weise 1991;Klevansky 1992;Volkov 1993;Hatsuda & Kunihiro 
1994; Buballa 2005). Recently, the trapped neutrino effect on the structure of neutron star 
was investigated in this model(Menezes,Providencia,& Melrose 2006; Ruster et al. 2006). 
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We describe the equation of state of the start with the two flavor NJL model defined by the 
Lagrangian density(Vogl & Weise 1991;Klevansky 1992;Volkov 1993;Hatsuda & Kunihiro 
1994;Buballa 2005) 

C = q (i^df, - m + #70) q + G [(qq) 2 + (qi^rq) 2 ] , (1) 

where m is the current quark mass, G the effective coupling constant, and fx the chemical 
potential matrix in flavor space, p, = diag (/i u , fid) = diag (/ig/3 — 2/x e /3, /Ab/3 + At e /3) with 
Hb and ji e being the baryon and electron chemical potential, respectively. 

The essential quantity characterizing a system in a grand canonical ensemble can be 
taken to be the thermodynamical potential Q. In mean field approximation, it can be 
expressed in terms of the effective quarks(Zhuang,Hufner,& Klevansky 1994), 

(m q - m ) 2 f d 3 p 



-e/^K (2) 



+Tln [(l + e -( E «+*0/ r ) (1 + e -(^-^)/r)] 
+Tln [(1 + e -(^+^)/ T ) (1 + e-^-^/T^ 

where m q = m — 2G (qq) is the constituent quark mass, (qq) the order parameter of the 
chiral phase transition, and E g = ^/va 2 q + p 2 the effective quark energy. To ensure the 
electric charge neutrality, we should include electrons in the system. Taking into account 
their contribution to the thermodynamic potential, 

- 2T /(05" 1 [( 1+e ^" ,/T » 
x (1 + e-< E '+<*>/ T ) 1 (3) 

with the electron energy E e = >/ml + p 2 and electron mass m e , the total thermodynamic 
potential of the system is 

n = n q + n e . (4) 

The physical order parameter (qq) or the effective quark mass m q should correspond to 
the minimum of the thermodynamic potential, 

d 

—n(T,n B ,n e ,m q ) = 0. (5) 

With the known thermodynamical potential, the pressure P and the energy density e 
of the system are related to f2 by 

p = -n, 

on on dn 
6 = - p - T &r-^Wu~* d Wa (6) 
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where we have considered the electric charge neutrality condition, namely that we choose the 
electron chemical potential /i e such that the system has zero net electric charge(Huang,Zhuang,& 
Chao 2003), 

d 

n e = -—Q(T,fi B ,fi e ,m q ) = 0. (7) 

The coupled set of gap equation and neutrality equation determines self-consistently the 
quark mass m q and the electron chemical potential \x e as functions of temperature T and 
baryon chemical potential \xb- Since we focus on compact stars, the temperature effect can 
be safely neglected. 

To simplify the calculation, we consider in the following chiral limit of the NJL model 
with mo = 0. There are only two parameters in the model, the coupling constant G and the 
hard momentum cutoff A to regulate the model. By fitting the pion decay constant and chiral 
condensate in the vacuum, they are fixed to be G = 5.02 GeV -2 and A = 653 MeV(Zhuang et 
al.1994). With these values, the effective quark mass m q and electron chemical potential /x e 
are shown in Fig {1] as functions of \i = /ie/3 at zero temperature. Since m q can be considered 
as the order parameter of chiral phase transition, its jump from 294 MeV to zero at \i c = 332 
MeV means a first order chiral phase transition. The system is in the chiral symmetry 
breaking phase at \x < fi c and the symmetry restoration phase at /J > fi c . The corresponding 
critical energy density at the transition is e c = 0.3 GeV/fm 3 . Due to the electric charge 
neutrality constraint, /i e jumps up at the transition point and is approximately proportional 
to fi in the symmetric phase. 
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Fig. 1. — The effective quark mass m q (solid line) and electron chemical potential /i e (dashed 
line) as functions of /x = /xb/3 at zero temperature in chiral limit of the NJL model. 

With the above equation of state with chiral phase transition, the structural of a non- 
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rotating compact star can be obtained by integrating the TOV equations(Glendenning 2000), 

dP_ (e + P) (M + 4?rr 3 P) 

~dV ~ N r (r - 2G N M) ' 



where = 6.707 x 10~ 39 GeV -2 is the universal constant of gravitation and M(r) the 
gravitational mass enclosed within the radius r. The chiral properties of the star is reflected 
in the pressure P(r) and energy density e(r). Giving a central pressure P at r = with /j, > /j, c 
which ensures that the center of the star is in the chiral restoration phase, one obtains the 
r-dependence of P, and the radius R of the star is determined by the condition P(R) = 0. If 
the star is assumed to be completely in the chiral restoration phase, one needs a bag constant 
to balance the nonzero pressure on the surface of the star. Taking the bag constant B = 75 
MeV/fm 3 , the maximum radius of the star is R = 7.1 km corresponding to the initial energy 
density e(0) = 1.5 GeV/fm 3 and final energy density e(R) = 0.6 GeV/fm 3 > e c . In our case 
with chiral phase transition inside the star, the core of the star is in the chiral restoration 
state and it is surrounded by the chiral breaking state. Taking into account the fact from the 
lattice simulations that the deconfinement and chiral restoration happen at the same critical 
point, the core can be regarded as a quark matter and the surrounding a hadron matter. 
For the initial energy density e(0) = 0.32 GeV/fm 3 , the chiral phase transition occurs at the 
critical radius R c = 5.2 km. 



3. Neutrino Oscillations 

The neutrino oscillation in matter is very different from the oscillation in vacuum 
(Caldwell,Fuller,& Qian 2000; Wolfenstein 1978;Mikheyev & Smirnov 1985;Rosen & Gelb 
1986;Bilenky & Petcov 1987; Gelb,Kwong,& Rosen 1997; Barget, Phillips, & Whisnant 1986; 
Kuo & Pantaleone 1989; Qian & Fuller 1995; Savage, Malaney,& Fuller 1991; Notzold & Raf- 
felt 1988; Botella et al.1987). We discuss in this paper only the v e — > v T oscillation in a 
compact star with chiral phase transition. In this case, the neutrino propagation in the star 
can be expressed as (Kuo & Pantaleone 1989) 

A - A cos 26 A sin 26 \ n / v e 
A sin 26 -A + A cos 26 ) \ \ v T 

where S = m% + m\ and A = m\ — m\ are related to the neutrino mass eigenvalues mi 
and m 2 in vacuum, and 6 is the mixing angle between the neutrino mass eigenvalue and the 
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weak interaction eigenvalue. Since £ contributes only a global phase factor to the neutrino 
eigenstate and does not affect the neutrino conversion (Fukugit a & Yanagida 2003), we omit 
it in the following. The matter dependence of the conversion is reflected in the factor A 
defined as 

A = 2V2G F EN e , (10) 

where E is the neutrino energy, Gf = 1.166 x 10 -5 GeV -2 the Fermi constant, and N e the 
free electron density 

d 3 p 1 



N ^ = 2 J (2^ i + e (^(r))/r - ( U > 

the r— dependence of the electron chemical potential is given by the TOV equation together 
with the equation of state. Suppose only electron neutrinos are created at the center of 
the star, the combination of the propagation equation (jSJ), the TOV equation (JSJ) and the 
chiral equation of state (j5J) controls self-consistently the space evolution of the transition 
from electron to tau neutrinos. 

To find the most suitable electron density N e where there is remarkable neutrino oscil- 
lation, we first consider the case with space-independent fi e . In this case, the propagation 
equation (jUJ) can be analytically solved, and we obtain the explicit surviving probability 

P e (r) = 1 - -sin 2 29 rn (l -cos^-M , (12) 



2 V A 

where sin 2 29 m and A m depend on fi e via 

(A sin 29) ' 



sin 2 26 r , 



{-A + A cos 29) 2 + (A sin 26 f ' 



X m = (13) 
v /(Asin2^) 2 + (A- Acos2#) 2 



In vacuum with N e = 0, the surviving probability is reduced to 

P e (r) = l-lsin 2 2^1-cos^. (14) 

It is easy to see that the biggest transition in matter happens at sin 2 29 m = 1 which 
corresponds to the matter environment with resonant electron density 

A cos 29 . . 

N " ^ 2^1' <15) 

where A, 9 and Gp are parameters fixed in the vacuum, the only adjustable parameter is 
the neutrino energy E. For low energy neutrinos, the resonant electron density N er , or 



- 7- 



equivalently, the resonant electron chemical potential /i er is large, and the matter induced 
oscillation is remarkably different from the oscillation in vacuum. With increasing neutrino 
energy, N er approaches to zero, and the difference in the oscillation between in matter and 
in vacuum disappears gradually. 

While the resonant electron chemical potential \i eT is independent of the chiral thermody- 
namics and the star structure, the matter effect to the neutrino oscillation is sensitive to the 
initial condition of the TOV equation, namely the equation of state at r = 0. For a smooth 
r— dependence of // e (r) around \i er characterized by the TOV equation, the r— integrated 
matter contribution to the neutrino oscillation is visible, while a sharp \i e {f) around fi er will 
make the matter effect to the neutrino oscillation difficult to be seen. 

Taking the vacuum parameters A = 2.6 x 1(T 3 eV 2 and sin 2 6 = 0.03(Maltoni 2004;Apol- 
lonio 1999), the effective neutrino mixing angle in matter is shown in FigfS]as a function of 
radius r for low energy neutrinos with E = 0.1 MeV. For the initial energy density e(0) = 0.32 
GeV/fm 3 , the peak sin 2 2# m = 1 is located at the resonant radius R r ~ 35 km > R c in the 
chiral symmetry breaking phase. This means that a remarkable matter induced neutrino os- 
cillation can not happen in a compact star without chiral phase transition, since the electron 
density N e in chiral restoration phase is too large and far away from the resonant electron 
density N er . 

To see the initial energy density dependence of the neutrino oscillations, we plot the 
electron chemical potential fi e around function of r in Fig[3]for different initial energy 

densities. The slope is very sensitive to the value of e(0). For e(0) = 0.32 GeV/fm 3 , \i e varies 
smoothly in the neighborhood of /j, er , and therefore, the space integrated matter contribution 
to the neutrino oscillation can be seen clearly. For e(0) = 0.39 GeV/fm 3 , \x e changes very 
fast around fi er , and the matter contribution occurs only in a very thin spherical shell of the 
star. In this case, the r— integrated matter effect is negligible. 

In FigJH we show the surviving probability (solid lines) of the electron neutrinos and 
their transition probability to tau neutrinos (dashed line) as functions of r. From our as- 
sumption of no tau neutrinos at the center of the star, the surviving probability P e (r) starts 
propagation with the initial value P e (0) = 1. In the core with chiral symmetry restoration, 
there is almost no neutrino conversion due to the too dense electrons. At the critical radius 
R c where the phase transition from chiral restoration to chiral breaking happens, the electron 
chemical potential and in turn the electron density drops down suddenly to a small value, 
and the visible conversion starts gradually. The biggest change occurs around the resonant 
radius R r . For low energy neutrinos with E = 0.1 MeV, almost all the electron neutrinos 
are converted into tau neutrinos around R r , and this conversion is kept in the further prop- 
agation. Clearly, the conversion probability is P T (r) — 1 — P e (r). With increasing neutrino 
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Fig. 2. — The effective neutrino mixing angle sin 2 29 m as a function of radius r of the star 
with initial energy density e(0) = 0.32 GeV/fm 3 for low energy neutrinos with energy E — 0.1 
MeV. 




r-R (km) 

Fig. 3. — The electron chemical potential /i e scaled by its resonant value ji er as a function of 
r — R r with R r being the resonant radius where the biggest neutrino transition occurs. The 
solid and dashed lines correspond, respectively, to the initial energy density e(0) = 0.32 and 
0.39 GeV/fm 3 . 

energy, the conversion is suppressed by the too large resonant radius in the chiral breaking 
state, as we explained above. For E ~ 0.3 MeV, the maximum change from electron to 
tau neutrino is reduced to about 50%, and for E > 0.8 MeV, there is almost no difference 
between the neutrino oscillations in matter and in vacuum. 
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Fig. 4. — The surviving probability P e of electron neutrinos (solid line) and their conversion 
probability P T to r neutrinos (dashed line) as functions of r for initial energy density e(0) = 
0.32 GeV/fm 3 and neutrino energy E — 0.1 MeV. 



We have examined the v e — * v T neutrino oscillation induced by chiral phase transition 
in a compact star. Due to the constraint of electric charge neutrality, the chiral properties 
controlled by baryon chemical potential fiB and the neutrino conversion governed by electron 
chemical potential /x e are related to each other. The star contains a dense core with chiral 
symmetry restoration rounded by the chiral breaking phase. When neutrinos propagate in 
the radius direction, the transition from chiral restoration to chiral breaking inside the star 
leads to significantly oscillations for low energy neutrinos. When the neutrino energy is low 
enough, the conversion probability is nearly 100%. Therefore, it is possible to detect the 
inner structure of compact stars in terms of low energy neutrino oscillations. 

We emphasize that the remarkable oscillation is resulted from the first-order chiral phase 
transition. For a second-order phase transition, due to the lack of the sudden drop of the 
electron chemical potential, the resonant radius R r will be several hundreds or even several 
thousands kilometers which are certainly beyond the reasonable radius of a compact star. 
In our calculation, R r is about 35 kilometers for initial energy density e(0) = 0.32 GeV/fm 3 
and is reduced to about 13 kilometers for e(0) = 0.39 GeV/fm 3 . Since the first-order chiral 
phase transition is a general consequence of lattice simulations and model calculations at 
high baryon density, our results obtained in the frame of NJL model is still qualitatively 
correct in general case. 

When other phase transitions at high baryon density are taken into account, what is their 
effect on the neutrino oscillations? Since the other transitions such as color superconductivity 
occur at extremely high density where the chiral symmetry is already restored, the electron 



4. Summary and Discussion 
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chemical potential \i e in these new phases is much higher than the resonant value /i er , and 
therefore, these new phases will not lead to remarkable changes in the neutrino oscillations. 
Similarly, the nonzero current quark mass m will change the behavior of [i e in the chiral 
restoration phase considerably, but it changes the neutrino conversion slightly. The neutrino 
oscillation in compact star is characterized by the first-order chiral phase transition. 

Acknowledgments: The work was supported in part by the grants NSFC10575058, 
10425810, 10435080 and SRFDP20040003103. 
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